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We show that the possibility to address and image single sites of an optical lattice, now an 
experimental reality, allows to measure the frequency-resolved local particle and hole spectra of a 
wide variety of one- and two-dimensional systems of lattice-confined strongly correlated ultracold 
atoms. Combining perturbation theory and time-dependent DMRG, we validate this scheme of 
lattice-assisted spectroscopy (LAS) on several example systems, such as the ID superfluid and Mott 
insulator, with and without a parabolic trap, and finally on edge states of the bosonic Su-Schrieffer- 
Heeger model. We highlight important extensions of our basic scheme to obtain an even wider 
variety of interesting and important frequency resolved spectra. 

PACS numbers: 


Probing time- or frequency-dependent response func¬ 
tions of strongly correlated systems forms the practical 
foundation for much of the current study of condensed 
matter. Among these, single particle observables have a 
special role, often constituting the basic measurements 
for determining a materials properties. Even individ¬ 
ual techniques in this field (STM, ARPES) form entire 
subbranches of experimental and associated theoretical 
physics of their own mm- To have such fundamentally 
important techniques also available for ultracold atomic 
gases has been the subject of intense research activity. 
For fermionic atoms, initial work focussed on gases in 
parabolic traps, first on obtaining the gap for paired 
fermions [3] by frequency-resolved measurements, then 
expanded to obtain the fully momentum-resolved spec¬ 
tral function via Bragg spectroscopy HE]. Performing 
this on systems that only experience an overall confin¬ 
ing potential but no additional optical lattice is due to a 
simple reason: typically, numbers in lattice-confined sys¬ 
tems are too small to obtain a significant signal. Also, 
such global measurements in presence of a trap come at 
a price, as the obtained Greens function may have their 
momentum-dependencies modified [6]. 

Yet, the presence of an optical lattice is key for the 
many proposed uses of ultracold atomic gases to serve as 
quantum simulators of strong correlation physics [7] , and 
practical ways for measuring time/frequency-resolved ob¬ 
servables in such systems have to be developed. Apply¬ 
ing Bragg spectroscopy, as has been done for ID bosonic 
Mott-insulators (MI), sidesteps any potential problems 
with insufficient atom numbers by performing the mea¬ 
surements on 2D arrays of decoupled ID chains, to obtain 
the spectrum of density-density [S] or single-particle ex¬ 
citations |9]. Performing Bragg spectroscopy in parallel 
in this way would however not be practical for 2D sys¬ 
tems, and further presupposes that the physics of each 


ID chain is not too affected by the inevitable density 
variations across the 2D array of decoupled chains. An 
alternative probing proposal, based on using a confined 
ion as a ’tip ’ m is still in practical infancy (though the 
setup has been partially realised m) 

Critically, the Bragg approach will fail for the latest 
generation of quantum gas microscope experiments [H- 
EI]. These may offer both single-site resolution and ad¬ 
dressing, but only across one or two parallel 2D planes, 
for what will essentially be systems of a few thousand 
atoms at most. Yet, the future study of strongly corre¬ 
lated atomic gases is moving very much in this direction, 
with multiple groups pushing to build up quantum gas 
microscopes. On the other hand, single-site addressing 
offers a very different alternative approach to measure 
frequency-resolved single particle excitations. Some of 
us [18] co-authored a recent proposal to leverage these 
new experimental capabilities to measure arbitrary time- 
dependent spin-spin correlators for any Hamiltonian im¬ 
plemented in an optical lattice that maps to an effective 
spin Hamiltonian. 

In this work, we show how to use single-site resolution 
and addressing to obtain time-dependent local particle 
and hole Greens functions for any Hamiltonian that can 
be prepared in a quantum gas microscope setup. We call 
this scheme lattice-assisted spectroscopy (LAS) in the fol¬ 
lowing. Using time-dependent density matrix renormal¬ 
ization group (t-DMRG) [T9| simulations we quantita¬ 
tively validate the LAS-scheme on both superfluid (SF) 
and MI systems, as well as on topologically protected 
edge states of the bosonic Su-Schrieffer-Heeger (SSH) 
model. We conclude with highlighting some of the most 
relevant of the many possible extensions of LAS. 

We consider the situation shown in Fig. [T] a: the sys¬ 
tem, assumed to be in the equilibrium state of some lat¬ 
tice Hamiltonian H (for which we want to measure lo- 
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FIG. 1: (Color online) a) Schematic overview of spectroscopic 
scheme: an extra outcoupling ‘site’ (OS) (physical or inter¬ 
nal state) is tunnel-coupled with amplitude J±{t) to lattice 
site X. The time-evolving average (hos(l)) on OS is obtained 
from repeated occupation-number imaging, b) Local den¬ 
sity for system in test case (A), 15 bosons in a ID parabolic 
trap. The MI region is in the centre, and SF regions occupy 
the wings. Red arrows indicate sites where system is locally 
probed, c) Symmetrized discrete absorption spectrum in cen¬ 
tre of MI region, extracted from the full time-evolving average 
(hos(i))(i“*^^ for J_l/J = 0.05 (black O), = 0.1 (red □), = 0.2 
(green o), and compared with {—u) obtained di¬ 

rectly from DMRG (blue solid line), d) (hos(l))(5“*^^^ over 
time at first resonance uj = 6.62J for J_i/J = 0.05 (black 
solid line), = 0.1 (red dotted), = 0.2 (green dashed). Rabi- 
oscillations can be clearly observed at experimentally acces¬ 
sible timescales, e) (fios(^))))“*^^ over time below first res¬ 
onance, at ai = 6J, vales of Jx/J match d). Off-resonance, 
oscillations lack all hallmarks of Rabi dynamics, exhibiting 
small amplitudes scaling with Jx/T, while frequency does not 
depend on it. f) Equivalent figure to c) for site 8 in the SF 
region. 


cal Greens functions), as well as an ‘extra’ site, denoted 
as the outcoupling ‘site’ (OS) in the following. The OS 
could denote another internal state of the atoms into 
which they are weakly coupled only on site a;, exploit¬ 
ing site-dependent addressing. Experiments have already 
demonstrated this capability m - thus all required tech¬ 
niques to implement LAS already exist. Alternatively, 
the OS could be another physical site initially separate 
from the system, with controllable tunnel-coupling to 
system-site x. Such a setup should be readily realisable 
with existing techniques - digital micromirror arrays inte¬ 
grated into quantum gas microscopes have already been 
used to write a large variety of spatially varying poten¬ 
tials resolved on the single-site scale [50]. With this, an 


OS could be created by allowing weak tunnelling from a 
system site only onto a single site in an empty parallel 
lattice, by imprinting a large energy-offset on every other 
site of the empty lattice. 

In either case, we assume an oscillating tunnel-coupling 
between atoms on site x and atoms in the OS. Its time- 
dependency may either be of form J±{t) = J±cos{uit) 
or J±{t) = J±sm{(jjt), switched on at time t = 0. Now, 
the resultant time-evolution of the average occupation on 
the OS will depend on the initial state of the OS. In the 
following, we focus on the most practically relevant cases, 
where there is initially either no atom on the OS, 


(nos(i))L“‘=° 

or exactly one: 
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In eq. (11): the first term is only present for bosons and 
absent for fermions. We also introduced the shorthand 
{nos{t))uj = {nos{t))cosu, + {nos{t))sinui- We assume that 
site OS has internal dynamics described by the Hamilto¬ 
nian hos = hosiji), diagonal in the occupation number 
n, with Shos{n,m) := hosin) — /los(w)- It is also as¬ 
sumed that only a single species of atoms may tunnel 
between sites x and OS, but LAS-formulas for multi- 
component systems are straightforward to derive. 

Here, the time-evolving averages are computed by con¬ 
sidering tunnelling as a quadratic perturbation, which 
links them to the local emission [absorption] spectrum of 
the system: 

p-PErr, 

n,m 

(3) 

integrated over with the kernel 


^t(w,cc') = 2^ 

(T —± 


(1 — cos((a;' -b auj)t)) 
{uj' + aujY 


(4) 


where |n), En, Z denote the exact eigenvectors, eigenen- 
ergies and partition function of the system Hamiltonian 
H. This kernel has the known property JCt{uJ,uj')/t 
(5(a; —CLi')-|-(5(a;-|-w'). If the probing time is not effectively 
infinite w.r.t. all other time-scales and the spectrum is 
not dense however, the implications of the kernel proper¬ 
ties for eqs. (0, (ll) fall broadly into one of two possible 
regimes: the excited states in'I^^^'°\u}') appear dense to 
lCt{uJ,Lj') at probing frequency uj and time t, or alterna¬ 
tively only one state lies close or at frequency uj inside the 
kernel envelope. In the first regime, the occupation on OS 
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will change linearly, {nos{t))ui {Ix{oj) +Ix{-oj))t 

- thus, the symmetrized can be extracted from 

a linear fit. In the second regime, at short times 
when the kernel is still very broad, {nos(t))uj will 
also behave linearly, as other eigenstates may still con¬ 
tribute. For longer times however, when ICt{ui,uj') nar¬ 
rows enough to filter out the discrete level, {nos{'t))oj 

W {En,rn:E^-E^=±u\{n\0\rn)\^s^^y^-, here 6 = 
aa;[aj.]. At even longer times, as the population of atoms 
on the OS begins saturating to the physically permitted 
level, we expect (and observe; see below) {nos{t))uj to be¬ 
gin deviating from the perturbative prediction strongly. 
On this largest timescale in the second regime, {nos{t))ui 
performs what are essentially Rabi oscillations, with fre¬ 
quency equal to 

at small enough J±, while the perturbative expression 
will naturally keep behaving monotonically. As the range 
of validity in time t and/or Jj_ is a key issue with any for¬ 
mulas from time-dependent perturbation theory we use t- 
DMRG in the following to validate LAS, and specifically 
eqs. Q and ([^, by comparing them to the full evolution 
of {nosit))u:- (where is determined independently, 

here also via DMRG). 

When experimentally applying LAS, both of the above 
regimes can be relevant. For one, cold atom experiments 
offering single-site addressability currently involve a few 
tens of sites and particles. At sufficient frequency res¬ 
olution (i.e. sufficiently long probing time) their bulk 
spectra will reveal their discrete nature. The other pos¬ 
sibility (applicable even when future experiments are so 
large that bulk spectra always appear dense), are local¬ 
ized gap-protected states, such as (topologically) pro¬ 
tected edge states or states around some impurity. As 
a model Hamiltonian to study the validity of eqs. Q,© 
for all these possibilities we use a generalised ID Bose- 
Hubbard model: 

L-l 

H = + h.c.) 

X — 1 

u ^ ^ 

+ — ^ fixifix - 1) + - {L + l)/2f fix, (5) 

X—1 X—1 

where aj, creates a bosonic atom on site x of an L-site 
chain with lattice constant a (with open boundary con¬ 
ditions) , tunneling between sites x and a: -|- 1 takes place 
with amplitude J -I- (—1)“(5J, V/ denotes the strength of 
any parabolic trapping potential that might be applied, 
and U > 0 denotes the strength of the onsite repulsion 
between atoms. 

We treated three test cases to study and demonstrate 
the capabilities of LAS. Test (A): Bosonic atoms in a 
parabolic trap form a so-called ’wedding cake’ when re¬ 
pulsive interactions are sufficiently strong m (c.f. Fig.§. 
For (5J = 0, Ip > 0, C/ J MI areas will alternate with 


SF transition regions from the traps centre outward. The 
LAS-scheme allows probing of the gapped or non-gapped 
nature of these respective regions directly. We choose a 
system of direct relevance to current experimental capa¬ 
bilities, 15 bosons in a trap with Vp/J = 0.06, U = 10 J. 
Its absorption spectrum can be obtained with great 
precision from quadratic fits to (hos(i))i/'*^^^ at short 
times together with observation of the Rabi oscillations 
at longer times (still occurring on experimentally acces¬ 
sible timescales), as outlined above, and summarised in 
Figs. c-f. This directly reveals the size of the gap in 
the MI centre of the ’cake’, and its disappearance in the 
wings. When the OS is coupled to a site in the MI and 
SF region respectively, we use the possible control over 
h{n) = Vpa'^{x— (L-I-1)/2)^ -|- [/os/2n(n— 1) to shift the 
contribution of to the emission spectrum, eq. (§ to 
be effectively zero, where Uqs denotes the atom-atom in¬ 
teraction on OS. In the MI centre, this is done by having 
Uos = 0, while in the SF wing this requires Uqs = U. We 
stress that even for a quite nonperturbative J± = 0.2J 
only very modest signal broadening is observed. 

For Test (B), we consider somewhat larger systems 
than previously, such as might be experimentally feasible 
in the near future. Here we want to show how LAS can 
be applied to obtain spectra that effectively approximate 
those of infinite systems, with a frequency resolution set 
by the finite size. For this, we focus on a Bose-Hubbard 
model with L = N = 51, SJ = Vp = h{n) = 0 in both the 
MI and SF regimes. We obtainin the emission spectrum 
I®™(w) from (hos(0)™*^°i from an OS coupling to the 
central site - note that we do not observe a gap in the 
MI regime as we are looking at the spectral function of a 
hole, and not of a particle (that would be For 

times t smaller than those at which the excitations trig¬ 
gered by the probing hit the systems open boundaries, we 
find excellent agreement with the finite-resolution spectra 
obtained from the convolution J du'dtlCt{u},uj')X^{u)'), 
which thus have a resolution set by the same time-scale 
(see Fig. [^. 

With Test (C), we highlight the power of LAS to 
probe the edge states present in topologically protected 
systems. For this, we study a system with L = 2N = 50, 
5J = 2/2>J, Vp = 0, h{n) = 0. The bosonic version of the 
SSH model [5T] has edge states near zero energy while the 
bulk is gapped (as would be expected of a model like ([^ 
with finite dimerisation, |(5J| > 0) as long as the particle- 
hole gap of the system is finite [32]. As shown in Fig. 
LAS applied to both edges and bulk allows accurate ob¬ 
servation of the distinct excitation spectra of both edges 
and bulk for larger U, and the gradual collapse of the 
distinct edge modes once U < 4( J -|- 6J). 

Our quantitative DMRG results also put to rest a po¬ 
tential key issue for LAS, the degree to which changes 
in the magnitude {nos(t))ui can be resolved in practice. 
With differences on the order of a few percent being read¬ 
ily detectable m. in every test case the largest J± value 
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FIG. 2: (Color online) Finite-resolution emission spectra in 
the thermodynamic limit for MI and SF phases, reconstructed 
from for J_i/J = 0.05 (black O), = 0.1 (red 

□), = 0.2 (green o), and compared with T^{lS) -|-I^™(—w) 
obtained directly from DMRG (blue x, blue solid line is guide 
to the eye) a) [/ = lOJ. h) U = 2J. 



FIG. 3: (Color online) Bulk and edge spectra of bosonic SSH 
model with L = 50. Markers denote spectra reconstructed 
from (nos(f))(((‘*^° for three different Jx/J: J±/J — 0.05 
(black O), J±/J = 0.1 (red □) and J±/J = 0.2 (green o). 
For comparison, we show I®™(a;)-|-I^™(—a;) obtained directly 
from DMRG (dashed blue lines for bulk at a; = 26 and solid 
purple lines for edge spectra, at a; = 1 and = 50, respectively). 
While U/{J + SJ) = 5, the distinct edge mode persists, while 
for lower U/{J -I- SJ) bulk and edge spectra start becoming 
indistinguishable, with overall spectral weight on the edges 
collapsing. 

we ran, J± = 0.2 J, always yielded clear signal wherever 
there was appreciable spectral weight, and in most cases 
J± = O.IJ gave useful data too. Even at Jj_ = 0.2J, 


spectral weight can be too small though - this is why we 
do not compare spectra from {nosit))ui against the low- 
magnitude parts ofIx{uj) in Fig.|^ As we show in Figs.[^ 
-[^the signal broadening at these Jj_ is modest - only for 
discrete excited state such as in the SSH model we find 
that amplitudes can be up to halved at J± = 0.2J when 
the spectral weight is low (yet still detectable). 

So far we have proposed and validated a general scheme 
to extract frequency-resolved local Greens functions for 
any system of interacting atoms confined to an op¬ 
tical lattice. Now we highlight several extensions to 
deal with nonlocal and higher-order correlations. First, 
if a single atom is created in a known superposition 
(ooQgi -b /3aos2)l®) between two OS (e.g. by having a 
single initially single-site localized atom delocalize ballis- 
tically, as in the so-called ’quantum hose-race’ [53]), it al¬ 
lows for the probing of non-local time-dependent particle 
Greens functions (a^(t)d|(0)). Specifically, when the sys¬ 
tem has inversion symmetry, the sites x and y (coupling 
respectively to OS 1, 2) are inversion-symmetric and 
/3 = ia and /iosi .2 = 0, we find (fiosi(0)<v-(^OS 2 (t))x = 
dt 2 J±{ti)J±{t 2 )lm[{ax-y{ti)al{t 2 ))]. In 
this way, the purely local contributions to (hosi, 2 )x 
are removed. Second, an atom initially prepared on 
an OS or two could be given a different internal hy- 
perfine state, making it distinguishable from the atoms 
in the system proper. Such schemes would provide a 
straightforward way to measure the local and nonlocal 
Greens function of a mobile distinguishable impurity in 
ID or 2D quantum systems, a subject which has at¬ 
tracted considerable theoretical [51H17I and experimen¬ 
tal interest recently [HIMSI]. Third, the cold atoms 
setup allows one to measure quantities that would be 
difficult to access in analogous condensed matter exper¬ 
iments, namely fourth-order time-dependent correlators 

(G) ^X3 ih) ^X2 {t2)axAti))- 
Measuring the time-evolving density-density correlator 
on OS 1 and 2, {hosi{t)noS 2 {t))uj (as in [ll|), coupling 
to lattice-sites x and y respectively with the same tunnel 
coupling, will yield 


{nosi{t)hos2{'t)) uj = 


1 




w I 


dt2 / dt[ / dt2J±{ti)J±{t2)J±{t2)J±{t'i) 

/o Jo Jo Jo 

^ \Cxyyx (tl, ^2 , ^2 j Cyxxy (G i ^2 , ^2 j ^ Cxyxy (Gjt 2 ,t 2 ,t]^) dl Cyxyx (G:t 2 ,t 2 ,t]^)] , ( 6 ) 


r 


where ± is the sign for bosons (fermions). As for eqs. 0. 
0. we assume only single-component tunneling to and 
from OS 1, 2. The detailed implications of these three 
extensions will be the subject of future work. 

In conclusion we have demonstrated a schemed to mea¬ 


sure correlation functions of ultracold atomic gases con¬ 
fined to optical lattices. This scheme only requires lo¬ 
cal addressability, which is available in today’s experi¬ 
ments. We have validated the scheme by comparing it 
to t-DMRG calculations and proposed extensions to deal 
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with non-local and higher order correlations. 

This work was supported in part by the Swiss NSF un¬ 
der division II, and the ARO-MURI grant (W9I1NF-14- 
1-0003). US acknowledges funding by the DFG through 
FOR801 and NIM. 


[1] A. Damascelli and Z.-X. Shen, Rev. Mod. Phys. 75, 473 
(2003). 

[2] O. y. Fischer, M. Kugler, I. Maggio-Aprile, C. Berthod, 
and C. Renner, Rev. Mod. Phys. 79, 353 (2007). 

[3] P. Torma and P. Zoller, Phys. Rev. Lett. 85, 487 (2000). 

[4] T.-L. Dao, A. Georges, J. Dalibard, C. Salomon, and 
I. Carusotto, Phys. Rev. Lett. 98, 240402 (2007). 

[5] J. T. Stewart, J. P. Gaebler, and D. S. Jin, Nature 454, 
744 (2008). 

[6] B. Frohlich, M. Feld, E. Vogt, M. Koschorreck, M. Kohl, 
C. Berthod, and T. Giamarchi, Phys. Rev. Lett. 109, 
130403 (2012). 

[71 1. Bloch, J. Dalibard, and W. Zwerger, Rev. Mod. Phys. 
80, 885 (2008). 

[8] T. Stoferle, H. Moritz, G. Schori, M. Kohl, and 
T. Esslinger, Phys. Rev. Lett. 92, 130403 (2004). 

[9] N. Fabbri, S. D. Huber, D. Clement, L. Fallani, C. Fort, 
M. Inguscio, and E. Altman, Phys. Rev. Lett. 109, 
055301 (2012). 

[10] C. Kollath, M. Kohl, and T. Giamarchi, Phys. Rev. A 
76, 063602 (2007). 

[11] C. Zipkes, S. Palzer, L. Ratschbacher, C. Sias, and 
M. Kohl, Phys. Rev. Lett. 105, 133201 (2010). 

[12] W. S. Bakr, J. 1. Gillen, A. Peng, S. Foiling, and 
M. Greiner, Nature 462, 74 (2009). 

[13] J. F. Sherson, C. Weitenberg, M. Endres, M. Cheneau, 
1. Bloch, and S. Kuhr, Nature 467, 68 (2010). 

[14] J. Simon, W. S. Bakr, R. Ma, M. E. Tai, P. M. Preiss, 
and M. Greiner, Nature 472, 307 (2011). 


[15] M. Cheneau, P. Barmettler, D. Poletti, M. Endres, 
P. Schauss, T. Fukuhara, C. Gross, 1. Bloch, C. Kollath, 
and S. Kuhr, Nature 481, 484 (2012). 

[16] T. Fukuhara, A. Kantian, M. Endres, M. Cheneau, 
P. Schaufi, S. Hild, D. Bellem, U. Schollwock, T. Gia¬ 
marchi, C. Gross, et al., Nat. Phys. 9, 235 (2013). 

[17] E. Haller, J. Hudson, A. Kelly, D. A. Cotta, B. Peaude- 
cerf, G. D. Bruce, and S. Kuhr, p. 6 (2015). 

[18] M. Knap, A. Kantian, T. Giamarchi, 1. Bloch, M. D. 
Lukin, and E. Dernier, Phys. Rev. Lett. Ill, 147205 
(2013). 

[19] U. Schollwock, Annals of Physics 326, 96 (2011). 

[20] P. P. J. Zupancic, Masters thesis, Ludwig-Maximilians- 
Universitat Munchen (2013). 

[21] W. P. Su, J. R. Schrieffer, and A. J. Heeger, Phys. Rev. 
Lett. 42, 1698 (1979). 

[22] F. Grusdt, M. Honing, and M. Fleischhauer, Phys. Rev. 
Lett. 110, 260405 (2013). 

[23] C. Weitenberg, M. Endres, J. F. Sherson, M. Cheneau, 
P. Schauss, T. Fukuhara, 1. Bloch, and S. Kuhr, Nature 
471, 319 (2011). 

[24] T. Kopp, A. Ruckenstein, and S. Schmitt-Rink, Phys. 
Rev. B 42, 6850 (1990). 

[25] A. Rosch and T. Kopp, Phys. Rev. Lett. 75, 1988 (1995). 

[26] M. B. Zvonarev, V. V. Cheianov, and T. Giamarchi, 
Phys. Rev. Lett. 99, 240404 (2007). 

[27] A. Kantian, U. Schollwock, and T. Giamarchi, Phys. Rev. 
Lett. 113, 070601 (2014). 

[28] S. Palzer, C. Zipkes, C. Sias, and M. Kohl, Phys. Rev. 
Lett. 103, 150601 (2009). 

[29] J. Catani, G. Lamporesi, D. Naik, M. Gring, M. Inguscio, 
F. Minardi, A. Kantian, and T. Giamarchi, Phys. Rev. 
A 85, 23623 (2012). 

[30] M. Koschorreck, D. Pertot, E. Vogt, B. Frohlich, M. Feld, 
and M. Kohl, Nature 485, 619 (2012). 

[31] C. Kohstall, M. Zaccanti, M. Jag, A. Trenkwalder, 
P. Massignan, G. M. Bruun, F. Schreck, and R. Grimm, 
Nature 485, 615 (2012). 



